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Abstract 

It is found that, in closed-Z-shell atoms, the exact local exchange potential v x (r) of the density functional 
theory (DFT) is very well represented, within the region of every atomic shell, by each of the suitably shifted 
potentials obtained with the non-local Fock exchange operator for the individual Hartree-Fock (HF) orbitals 
belonging to this shell. Consequently, the continuous piecewise function built of shell-specific exchange 
potentials, each defined as the weighted average of the shifted orbital exchange potentials corresponding to 
a given shell, yields another highly-accurate representation of v x (r). These newly revealed properties are 
not related to the well-known step-like shell structure in the response part of v x (r), but they result from 
specific relations satisfied by the HF orbital exchange potentials. These relations explain the outstanding 
proximity of the occupied Kohn-Sham and HF orbitals as well as the high quality of the Krieger-Li-Iafrate 
and localized HF (or, equivalently, common-energy-denominator) approximations to the DFT exchange 
potential v x (r). The constant shifts added to the HF orbital exchange potentials, to map them onto v x (r), 
are nearly equal to the differences between the energies of the corresponding KS and HF orbitals. It is 
discussed why these differences are positive and grow when the respective orbital energies become lower 
for inner orbitals. 

PACS numbers: 31.00.00, 31.15.E-, 31.15.xr 
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I. INTRODUCTION 



Representing the quantum state of a many-electron system in terms of one-electron orbitals is 
simple and theoretically attractive approach. Such description is realized in the Hartree-Fock (HF) 
method yj], as well as in the Kohn-Sham (KS) scheme of the density-functional theory (DFT) \2, 
dLca\- The latter is an efficient and robust tool which is now routinely applied in the calculations of 
electronic properties of molecules, even very large and complex, and condensed-matter structures. 
Though the KS scheme is formally accurate, the one-body KS potential contains the exchange- 
correlation (xc) potential v xc , whose exact dependence on the electron density remains unknown. 
It is usually treated within the local-density or generalized-gradient approximations (LDA, GGA), 
despite the well-known shortcomings of the LDA and GGA xc potentials (especially the self- 
interaction errors). Some of these deficiencies are removed when the exact form (in terms of the 
occupied KS orbitals) is used for the exchange part E x of the xc energy. The exact exchange 
potential v x is then found from E x by means of the integral equation resulting from the optimized- 
effective-potential (OEP) approach j^a)]^, [7, [si I9] or by using the recently developed method 



based on the differential equations for the orbital shifts [fiol [nil; another method based on the 



direct energy minimization with respect to the KS-OEP potential (expressed in a finite basis) I112h 



suffers from convergence problems [13] which are not fully resolved yet and they are still under 
study 11141 1 1511 . The exact potential v x is free from self-interaction and it has correct asymptotic 
dependence (— 1/r for finite systems) at large distances r from the system; thus, unlike the HF, 
LDA or GGA potentials, it produces correct unoccupied states. In the DFT, the approximation, 
in which the exchange is included exactly but the correlation energy and potential are neglected, 
is known as the exchange-only KS scheme — it is applied in the present investigation. The full 
potential v xc can also be found by means of the OEP approach when the DFT total energy includes, 
besides the exact E x , the correlation energy E c depending on all (occupied and unoccupied) KS 
orbitals and orbital energies [8]. This makes such computation tedious, to a level undesirable in the 
DFT, since it involves calculating E c with the quantum-chemistry methods, like the M0ller-Plesset 
many-body perturbation approach. 

Defined to yield the true electron density, the KS one-electron orbitals have no other direct 
physical meaning since they formally refer to a fictitious system of non-interacting electrons. 
However, it is a common practice to use these orbitals in calculations of various electronic proper- 
ties; in doing so the iV-electron ground-state wave function of the physical (interacting) system 
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is approximated with the single determinant built of the KS orbitals. This approximate approach 
is justified by (usually) sufficient accuracy of the calculated quantities, which is close to, or often 
better than, that of the HF results llql . It seems that the success of the DFT calculations would 
not be possible if the KS determinant, though being formally non-physical, was not close to the 
HF determinant which, outside the DFT, is routinely used to approximate the wave function ^ 
of the real system. Therefore, understanding this proximity is certainly very important for the 
fundamentals of the DFT 

V the whole KS and HF deter- 



minants 



Previous calculations [5(a) 



21, 



22, 



2311 have shown that, not on 



17|] and the corresponding electron densities [5(a) 



18, 



3Q], 



but also the individual 



occupied KS and HF orbitals, a <r( r ) and 0™(r), in atoms are so close to each other that they 
are virtually indistinguishable (here the orbitals, dependent on the electron position r and the spin 
a —{, f , are numbered with index a = 1, . . . , N a ; N± + iVj = N). This property is particularly re- 
markable for the exchange-only KS orbitals which differ so minutely from the HF orbitals that, for 



18] 



atoms, the OEP total energy is only several mhartrees higher than the HF energy _E H f PLIES. 
The outstanding proximity of the KS and HF orbitals is surprising in view of the obvious differ- 
ence between the exchange operators in the KS and HF one-electron hamiltonians (see below) and 
the fact that the corresponding KS and HF atomic orbital energies, e a „ and e™, differ substan- 
tially, up to several hartrees for core orbitals in atoms like Ar, Cu 1(1 §1 [except for the KS and 
HF energies of the highest-occupied molecular orbital (HOMO) which are almost identical]. This 
apparent contradiction has not yet been resolved; in Ref. it is suggested that the KS and HF 
determinants are close to each other "since the kinetic energy is much greater than the magnitude 
of the exchange energy". 

The present paper investigates the proximity of the KS and HF orbitals and it reveals that, in 
closed-Z-shell atoms, there exists a direct mapping between the HF orbital local exchange poten- 
tials v™ a ( r ) an d the DFT exact local exchange potential v xa (r). The former are specific to each 



HF orbital 0™( r ) an ^ 316 defined as 



(1) 



with the Fock exchange non-local operator v^ a (r) within the HF approximation that describes the 
interacting system. The DFT exchange potential v xa (r) is common for all orbitals relevant to the 
KS non-interacting o subsystem. This potential is found to be very well represented, within the 



region of each atomic shell, by the individual, suitably shifted potentials v"^(r 



v HF (r) 



obtained for the HF orbitals that belong to this shell; the constant shifts C aa are orbital-specific. 
As a result, for each shell, the weighted average of the potentials vf^ a (r) corresponding to the 
orbitals from this shell yields the shell- specific exchange potential that also represents v xa (r) with 
high accuracy within the shell region. The revealed mapping between v^ a (r) and v xa (r) is shown 
to have origins in the specific relations satisfied by the HF orbital exchange potentials. Thus, the 
proximity of the KS and HF orbitals is explained. Simultaneously, it becomes clear why, in atoms, 
the exact exchange potential v xa (r) (where r = |r|) has the characteristic structure of a piecewise 
function where each part spans over the region of an atomic shell and it has distinctively different 
slope dv xa (r) / dr in consecutive shells 1251] . 

The specific properties of v^ a (r) are also shown to be directly responsible for the high quality 
of the approximate representations of the exact exchange potential v xrT (r) that are obtained in the 
Krieger-Li-Iafrate(KLI) [5] and localized HF (LHF)[17] approximations, the latter of which is 
equivalent to the common-energy-denominator approximation (CEDA) |]26ll . The constant shifts 
C aa , needed to map the HF potentials v™ a (r) onto v xa (r), are shown to be nearly equal to e aa — e^. 
This leads to better understanding why, for each KS occupied orbital (other than the HOMO), its 
energy e aa is higher than the corresponding HF energy and the difference between these two 
energies is larger for the core orbitals than for the valence ones. Finally, it is shortly argued that 
the presently revealed properties of the KS and HF exchange potentials do not result from the 



we 



1-known step-like shell structure present in the response part v™ v (r) of the exchange potential 



121281. 



H. THEORY 



A. Hartree-Fock method and optimized-effective-potential approach 

The HF one-electron spin-orbitals 4>™( r ) 3X6 obtained by minimizing the mean value {^\H\^) 
where H is the Hamiltonian of the iV-electron interacting system and \& belongs to the subspace 
fijSf of normalized A r -electron wave functions that are single Slater determinants built of one- 
electron orbitals. Similar minimization is carried out in the exchange-only OEP method, but there 
is the additional constraint that for every trial determinant all N a constituent spin-orbitals </> a(T (r) 
satisfy the KS equation with some local KS potential v sa (r). The minimizing potential v SCT (r) = 
t> s ° EP (r), yields, after subtracting from it the external v ext (r) and electrostatic v es (r) terms, the exact 
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exchange potential v xcr (r) 
(fraa), so that we have 



,OEP, 



r) (corresponding to the density n a calculated from occupied 



(r) +f es (r) + v xa (r) . 



(2) 



It has to be stressed here that the proximity of the exchange-only KS and HF orbitals is not readily 
implied by the fact the two sets of orbitals result from the minimization of the same functional of 
energy, i.e., E\^>] = (ty\H\ty) where \I> G fijSf . Indeed, for a suitably chosen model Hamiltonian 
H, the corresponding HF orbitals 0^( r ) tnat minimize E[*ff] might not be well approximated by 
any set of one-electron (KS) orbitals acr (r) that come from a common local potential u s<T (r). Then, 
the latter condition, which is imposed on the orbitals (j) aa {r) in the OEP minimization, would be so 
restrictive that the obtained KS-OEP orbitals would differ significantly from the HF ones. Thus, 
it seems that it is the specific form of the physical Hamiltonian H (with Coulombic interactions) 
that actually makes the close representation of the HF orbitals with the KS ones possible. 

The exchange-only KS equation, satisfied by the corresponding (OEP) orbitals </> a(T (r) and their 
energies e acr , takes the form 

1. 



h s Jr) 



~2 V + v ^ 



r) + v es (r) 



r) (f) a a(r) = e 



(r) 



(3) 



(atomic units are used throughout) where we put v xa (r) 
electron density n tot (r) = n-f(r) + rij(r), which enters 

a ' Ut 
ar — 



v es [n tot ](r 

is the sum of the spin-projected densities 
In the HF equation 



j?F(r) in the OEP case. The total 



(4) 



a=l 



(5) 



+ v ext (r) + C(r) + vl(r) CW = « 



(6) 



satisfied by the orbitals </>^(r) and energies e™, the multiplicative local exchange potential v XfT (r), 
present in the KS equation ©, is replaced with the non-local Fock exchange integral operator 
^(r), built of {4>w}a=i' i ts action on a given HF orbital (f)™ yields yj] 



6=1 



dY , €W) 



(7) 



The electrostatic potential f^ F (r) = t> es [n™](r) is found for the HF total electron density n™(r) 
defined in a similar way as n to t(r). The KS and HF orbitals are ordered according to non- 
descending values of the corresponding orbital energies ( e acr < e a+ i )Cr and < e^ 1(T , 
a = 1, . . .)• Both the KS and HF equations need to be solved selfconsistently. Real KS and 
HF orbitals are used throughout this paper. 

Obviously, for each HF orbital <p™{r), the Fock exchange operator ^(r) present in the HF 
equation © can be formally replaced by v™ a {v), Eq. (OQ), however, this local exchange potential 
is orbital-dependent due the non-locality of ^(r), Eq. ©. Thus, also, the resulting total HF 
potential 

vfM = ««t(r) + C(r) + Crto (8) 

is different for each orbital <p™(r), unlike in the KS scheme where all electrons (of given spin 
a) are subject to the same total potential v sa (r), which includes the common exchange potential 
v xa (r). Dependence on a will be suppressed hereafter (unless otherwise stated). 

B. Orbital and energy shifts. Exact exchange potential 

The exact exchange potential v x = v x EP satisfies the OEP equation iRfloll 

N 

Mr) = 2^0 a (r)50 a (r) = O, Vr , (9) 

a=l 

which results from the OEP minimization and depends on v x through the orbital shifts (OS) 8(j) a (r). 
Each OS fulfills the equation fl fiol [llll 

h s (r)-e a \64 >a (r) = W a ± (r) (10) 

(where 4> a , e a are the solutions of Eq. ©) and it is subject to the constraint (0 a |50 a ) = 0. The 
equation (fTOl) includes the KS Hamiltonian h s , present in Eq. ©, and the term (defined using the 
sign convention of Refs. 

W^(r) = [vl(r) + D aa - v x {r)] a (r) . (11) 

where v x = v x [{(fib}] and 

D aa = (0aK ~ ^|^o) ■ ( l2 ) 

It should be noted that / dr0 a (r)W^-(r) = 0. 



The OS 5(f) a and the energy shift (ES) - the constant D aa give, within the perturbation theory 
(PT), the first-order approximations to the orbital and energy differences (shifts), — (0™ — (f> a ) 
and — (e™ — e ), respectively. Here, the orbitals 4>™ and the corresponding energies e™, are the 
solutions of the HF-like equation which is the same as Eq. © except for v x replaced by built 
of the KS orbitals {<pb}- The corresponding perturbation is then equal to 5h s = v* — t> x so that the 
first-order correction to e a is — 5e a = (0a|<5^ s |0a) = —D aa while the correction to a (r) is 

oo 

-6<p a (r) = ^2 c i«^( r )' ( 13 ) 

c ja = , (14) 

e j - e a 

-Dja = -{4>j\8K\4> a ) = {4>j\v x - vl\4>a) ■ (15) 

It satisfies Eq. (fTOl) and the constraint (0 a |<50 a ) = indeed. Obviously, the solutions 0™, 
are not identical to the selfconsistent HF orbitals (p™ and orbital energies which are obtained 
from Eq. ©. The latter HF quantities can also be found within the PT approach by calculating 
the differences A<fr a = (f)™ — 4> a , Ae a = e^ F — e a in the first-order approximation. In this case, 
the perturbation is given by Ah s = hup — h s [where hup is the HF Hamiltonian of Eq. ©] and 
it consists of three terms, Ah s = 5h s + Av es + A-0 F . The terms Av es = v es [n^] — v ef ,[n tot ] = 
v es [n™ — 7Ttot] (cf. Eq. ©) and Av* = i'l[{(j)a F }] — v*[{4>cl}] depend on A(j) a (of both spins for 
At> es ), linearly in the leading order, so that they have to be calculated selfconsistently even in the 
PT approach. But, if we substitute (— S(j) a ) for A(p a the difference — n tot becomes 5n^ + dn^ 
so that it vanishes due to the OEP equation ©. Then, we find Av es = and the perturbation Ah s 
becomes Sh s + A-0 F {{4> a }, {8<f>a}]- It can be further reduced to 5h s if the OS 5(p a are sufficiently 
small. This argument, although not strict, leads to the conclusion that the differences A(p a and Ae a 
are well represented by the orbital and energy shifts, — 8(j) a and — 5e a = —D aa , respectively, which 
are obtained with the perturbation 5h s . This conclusion is confirmed by the relations \\A(L — 
{-6<f> a )\\ < O.13||A0J (where ||0|| 2 = / dr\<P(r)\ 2 ) and |Ae a - (-<fe a )| < 0.003|Ae a | [U, 
established numerically for the Be and Ar atoms (see Tables U and HO); the above inequalities are 
obtained for <p a , e a , 5(p a calculated as in Ref. Jilt], and <p™ (expanded in the Slater-type-orbital 
basis), e^ F taken from Ref. J3p| . The representations of <p™ — <p a by — 5<p a and — e a by 
— 5e a will be used in further discussion. They can also be applied to construct a nearly accurate 
approximation of the exact exchange potential; the new method will be reported elsewhere soon 

bit. 
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The part of W a (r) = 5h s (r)(p a (r) parallel to the orbital <p a is 

WKr) = -^.(r) 



(16) 



and it sets the ES 5e a = D aa . The part 

W£(T) = W a (r)-W!(r), (17) 

perpendicular to <p a , sets the OS 5<p a {r), Eqs. (flOl) . (ITTb . Thus, the KS and HF orbitals, a (r), 
0™(r), can be close to each other, even if the orbital energies e tt , e^ F , differ significantly, provided 
the term W^"(r) is sufficiently small. Note that the orbitals remain unchanged when a (possibly 
orbital-dependent) constant is added to the Hamiltonian in the KS or HF equations. 

When the equation (flOl) (after multiplying it by a (r) and subsequent summin g oy er a = 
1, . . . , N) is combined with the OEP condition ©, the following expression UJl for the 

exact exchange potential is obtained 



-x OEP (r) = ^ KLI (r) + ^ x OS (r) (18) 
It contains the KLI-like potential [50 

K0a}, {D aa }] (r) = ^ x sl (r) + ^ x ES (r) (19) 
which consists of the Slater potential 

v x S1 W = ^TE^( r ) 5 x F W^(r) (20) 

and the ES term, linear in D aa , 

1 - 

^ ' o=l 

where n(r) = J2a=i ( r ) > these terms are defined with the OEP orbitals <p a (r) and constants 
D aa . The OS term present in Eq. (fT8~l) . linear in S<f> a (r), is 

1 - 

^x OS W = — -^[2e a a (r) - (V0 a (r)) • V]50 a (r) . (22) 

Since any physical potential is defined up to an arbitrary constant, it is usually chosen that the 
constant D NN = for the HOMO [10]; then the potential v° EP (r) goes to as — 1/r for r = 

8 



Irl — > oo (except for the directions that lie within symmetry planes in some molecules: in this 

nn 

special case the (— 1/r + const) dependence at large r is found; cf. Ref. U0U 32[). 

However, the use of Eq. (fl8l) for calculation of u° EP (r) still requires solving the equations 
(191 101) for 5(f) a (r) as well as determining the selfconsistent values of the constants D aa which 
depend on v x = t>° EP (r) through Eq. (PT21) . This solution is obtained in an iterative way in Ref. 
, while a non-iterative algorithm, where both sets {5(f) a } and {D aa } are found simultaneously, 
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is presented in Ref. [ll|. Let us note that the equations (|9lTT2l). (118H22I) can be used to determine 
the exact exchange potential t>° EP (r) not only in the exchange-only OEP approach, but also when 
the orbitals <j> a (r) are the solutions of the KS equation with the potential u s (r) that, besides v x (r), 



includes a correlation term v c (r) 



C. High-quality KLI and LHF (CEDA) approximations 

Since the OS 8<p a {r) are usually small, the term v® s (r), Eq. ([22]) . is a minor correction to 
f^P^r) in Eq. (TT8T) . Therefore, when we neglect t>x S (r) completely, the exact exchange potential 
■u° EP (r) is represented with high quality by the KLI-like term vf L1 (r), Eq. (fT9l . The original KLI 
approximation ||5|] 

*£" [{0a}] (r) = v^ 1 [{<f> a }, {D™}} (r) (23) 
is obtained (here for the KS-OEP orbitals a ) when the constants 

= (<Pa\vf U - V F x \<Pa) (24) 

are found selfconsistently, analogously as D aa in Eq. (fT2l) for v x = v° EF . Since, the equation (1241) 
remains satisfied when an arbitrary constant, but the same for all a, is added to each D^ 1 , one 
usually sets -D^v = which makes the potential t^^r) decay like —1/r for large r. 
The sum over j in Eq. (fT3l can be split into two terms, 

N 

= - E (25) 

b=l,b=ia 
oo 



which are the projections of the OS 5cf) a onto the subspaces of occupied (occ) and virtual (vir) 
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orbitals, respectively. Thus, the OS term v ° s (r), Eq. (|22l) . can be rewritten as follows 

JV-l N 



) = E E 2 ^0a(r)0 b (r) + ^ x os - r (r) (27) 



t>° s ir 

v 7 a=l b=a+l 



after the definition (1141) of Q a and relation D ba = D ab [cf. Eq. (1151) . i£ is Hermitian and real] are 
used; the term 

u OWr (r) = v OS [{0a}) {eJ) ^Jrj] (r) (2g) 

is found by substituting 5 (ft™ for m Eq. (1221) . When the OS <50 a are small, the corresponding 
projected parts 6(f)™ are even smaller since the general relation ||<50° cc || 2 + ||<50a ir || 2 = ll^0a|| 2 
holds. Then, another high-quality representation of v ° EP 

1 - 

£ X LHF [{04, {D ab }) (r) = t; x sl (r) + IW«(r)&(r) (29) 

^ ' a,b=l 

is obtained by setting <50™ r (r) = in the OS term t£^(r), Eqs. (|2"7ll28l) . This representation yields 



the well-known LHF (CEDA) approximation II 171 



26] 



^ LHF [{0j]=^ LHF [{0j,{^ L n] oo) 



(here defined for the set {4> a } of the KS-OEP orbitals) when the constants D]^ ¥ = (4> a \vf u — 
v ¥ \(j)b) , defined analogously as in Eq. (fT5l) . are found selfconsistently for (ab) ^ (NN); we 
also set D)$y = 0, as in the KLI case. Let us note that the condition 5(f)™ = is equivalent 
to the relation (p^. = acr + Ylb£ a c ba4>ba (valid in the first-order approximation) which, when 
satisfied for both spins a, implies that the HF determinant built of {0™} is identical to the KS 



17 



to derive 



determinant built of {<fi aa }. This (approximate) identity has been assumed in Ref. 
the LHF approximation. Obviously, both the KLI and LHF approximate exchange potentials can 
be defined for any set of (orthogonal, bound) orbitals {4> a }a=i- I n particular, it can be done for 
the orbitals that are selfconsistent solutions of the KS equation © where the potential v x is set to 
v™[{<p a }] or [{£,}]. 

The high quality of the KLI and LHF approximate potentials, when derived as presented above, 
clearly results from the proximity of the HF and KS-OEP occupied orbitals which is characterized 
by the small OS 5(p a . However, the OS terms v° s and t> x S vir which are neglected in the KLI and 
LHF (CEDA) approximations, respectively, are expressed through all OS 5(f) a (or their projected 
parts 5(j)™'). As a result, some information associated with the small magnitudes of the individual 
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OS 5 (pa may be lost in the resulting potentials and . In particular, the Slater term, Eq. 
(|2Q|) . present in these potentials, can be viewed the weighted average 

«?(■•) = X>»(r)^ (31) 

a=l ^ > 

of the KS orbital exchange potentials 

,.W = *P, (32 ) 



so that it cannot fully reflect the properties of the individual v xa (r). In the following discussion 
(Sec. InTl) for closed-/-subshell atoms, new properties of t> xa (r) are exposed on/j when the prox- 
imity of the HF and KS-OEP orbitals is considered separately for each orbital. 

D. Closed- /-subshell atoms: Fock exchange operator, orbital exchange potentials 

For a closed-Z-subshell atom, the non-local (integral) Fock exchange operator, acting on an 
atomic orbital a (r) = r~ l Xni(r)Yi m {6 ', ip) (a = nlm), yields 

£ x F (r)0 a (r) = r-V^r) Y lm (9, 0) (33) 

where Y[ m (9, (p) is the spherical harmonic, Hereafter, the orbitals are labeled with the principal, 
orbital, and magnetic quantum numbers, n, I, m; the symbols n occ and /max will denote, respec- 
tively, the largest number n and the maximum value of I for given n, within the set {occ} of the 
occupied orbitals {4> n i m } (hereafter, we refer to this set with the general label "occ"). It will be 
convenient to have a notation for the HOMO label: H = (n ImL) at n = n occ ; note that the HOMO 
belongs to the outmost occupied shell for the closed-/-shell atoms. The factor 

occ l+l' 

F x , nl (r) = J2 E 9{l,l',l")xMr)vvin'l',nl;r), (34a) 

n'l' l" = \l-l'\ 

is defined [1] (here with the occupied KS radial orbitals Xn'i'(r) ) through the functions 

f°° (r V" 

v l „(ril',nl;r) = - / dr' '^F+lXn' v (r')xm(r ') (34b) 



( r >) 
. 1 V l " \ 

where we denote g(l, I', I") = (2V + 1) (a special case of the 3j Wigner symbol), 

min(r, r'), r > = max(r, r'). In particular, the following non-zero coefficients g(0, 0, 0) = 1, 
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0(0, 1, 1) = 1, g(l, 1, 0) = 1, g(l, 0, 1) = 1/3, g(l, 1, 2) = 2/5 are needed to find the quantities 
F X ;ni( r ) for atoms with s and p orbitals (like Be, Ar); note that the step in the summation over I" 
in Eq. (|34a| ) is 2. Thus, the orbital exchange potential, Eq. (|32l . 



Vxa(r) = v x . nl (r) = F x . n i(r)/Xni(r) 



(35) 



is obtained; the corresponding HF quantities, denoted as v^F(n'l',nl;r), F x H ^(r), f^ a F (r) 
v™ nl (r), can be determined with the HF atomic radial orbitals Xnf( r )- The OS 



r 1 $Xni{r)Y lrn (6,(p) 



depends on the term 



through its radial part 



W a ± (r)=r- l W^(r)Y lm (9,v) 



WT'ir) = F x . nl (r) + [D nl , nl - v x (r)\ XnI (r) 



(36) 



(37) 



(38) 



entering the equation 



1 d 2 1(1 + 1) 



+ vJr) - t n l 



(39) 



2 dr 2 2r 2 

for Sxni(r) derived from Eq. (flOl) ; here e n i is the energy of the KS orbital <p a = 4> n i m . The KS 
potential v s (r), Eq. © contains the term f ex t(^) = —Z/r where Z is the atomic number, equal to 
iV for neutral atoms. 



III. NUMERICAL RESULTS AND DISCUSSION 



A. Proximity of KS and HF orbitals 



The proximity of individual HF and KS orbitals can be quantified with the norms ||50 a || which 
are found to be indeed very small, in comparison with ||0 O || = ||0„ F || = 1. Calculating the OS S(j) a 



with the method of Ref . 



11L we obtain \\S(f) a \\ < 0.007 for each occupied orbital in the Be and Ar 



atoms; see Table HI The partition 



\\5<t>nlm\\ 2 =J2 C n%nl, (40) 
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plotted for Ar in Fig. [fl shows that, among the KS bound orbitals (p n 'i m , the dominating contribu- 
tions (? n , ln[ to the nlm OS come from the n'lm orbitals with n' — n — 1 and/or n + 1, i.e., from 
the neighboring electronic shells; e.g., for S(j) 3s in the Ar atom, the largest terms c^, lnl are found 
for the n'l = 2s (occupied) and n'l = 4s (unoccupied) orbitals. But, there remains a large part 
of ||<50 n z m || 2 which cannot be attributed to higher unoccupied bound states 4> n 'i m since the corre- 
sponding (? nn . n i terms vanish rapidly with increasing n'\ see Fig. Q] This unaccounted part comes 
from continuum KS states {e n >i > 0). Let us also note that, for each OS 50 n ; m analyzed in Fig. [Q 
its projection 8<j>^ m , Eq. (|25l) . onto the occupied- state subspace has the squared norm smaller than 
\\\ 84>nim || 2 which means that the relation ||50°^|| < ||<50™mll holds for the Ar atom. 

The above results also confirm that the assumptions 5(p a = and 5(p™ = 0, which can be used 
to derive the KLI and LHF (CEDA) approximations, respectively (cf. Sec. Ill CI) , are very accurate 
but not exact. 

B. Exact exchange potential vs orbital exchange potentials 

The norms ||50 n z m || have such low values because the terms W^' rad (r) are sufficiently small 
for all r (the scale of this smallness will be discussed later on). This, combined with the relation 

v^ni (r) + D nl . nl = v x (r) + nl . \ } , (41) 

Xm{r) 

found with Eqs. (|38l) and (1351 ), implies that each shifted orbital exchange potential (calculated 
from the KS-OEP orbitals) 

v x -ni(r) = v x - n i(r) + D n i- n i (42) 

is very close to the exact exchange potential v x (r) = v® EP (r) within the r-interval 
( r n-i,n, r n,n+i) = S n where the denominators in the right-hand side of Eq. (I4TI) . i.e., the or- 
bitals Xni(r) from the n-th atomic shell (K, L, M, . . .), have largest magnitudes. The shell border 
points r n>n+ x for n = 0,1,..., n occ — 1 (the respective HF points r™ n+1 , defined precisely below, 
can be used) are near the positions r™ n where the radial electron density p(r) has local minima. 
In large parts of the shells S n <, n' < n, where the orbital Xni( r ) entering the denominator in Eq. 
(I4TT) has sizeable magnitude (though at least a few times smaller than in the shell S n ) the potentials 
v x - n i(r) are also close to v° EP (r) (but not so tightly as for r G S n ). The proximity of the poten- 
tials is evident in Figs. El El S] for the Be, Ar, and Zn atoms, respectively; it also holds for other 
closed-/-subshell atoms. It is disturbed in the vicinity of the nodes of Xni(r), where the potential 
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v x;n i( r ) diverges while the term W nl ' ra (r) is finite and small. The potential v X]n i(r) also differs sig- 
nificantly from r OEP (r) within the occupied shells S n ", n" > n, where both the functions Xni(r), 
W^ ;rad (r) decay exponentially. 

In the asymptotic region (spanning outside the occupied shells, i.e., for r > r n >n+1 , n = 
riocc) the exact exchange potential r> OEP (r) lies very close only to the HOMO exchange potential 
v X ;H( r ) = v x,H ( r ) (where D H H = 0) which has the correct — 1/r dependence for large r resulting 
from Eqs. (1341 ), (1351) ; see Fig. [5] Indeed, the potential v x - n i(r) for nl ^ H includes, besides 
the self-interaction term Vo(nl, nl; r), equal to —1/r for large r, also, at least one non-zero term 
proportional to Xh{^)^i" (H, nl;r)/xru( r ) with/" ^ 0;cf. Eqs. (|34l) . ((3~5l) . The latter term diverges 
for r — > oo since the factor vi»(H, nl; r) tends to a constant while each KS radial orbital Xni( r ) 
decays like r 1 ^ nl e~ l3nir where (3 n i = y/—2e n i (cf. Ref. Q); this is true also for nl = H. The Be 
atom, with the Is and 2s orbitals only, is the only exception here since, in this case, both potentials 
v x -u( r ) an d v x -2 S (r) decay as —1/r for large r. Indeed, with Eqs. (|34l ), (|35l ) we find the following 
expressions 

v x -u{r) = Vq(1s, Is; r) + X2s ^\ vo{2s, Is; r) , (43a) 

v x -,2s{r) = v {2s, 2s; r) + ^4^^ (ls, 2s; r) , (43b) 

valid for the Be atom. Due the orthogonality of the Is and 2s orbitals, the function r (ls, 2s; r) = 
v (2s, Is; r), Eq. (|34bl) . is equal to J r °° dr'(l/r — I / r')xi s ( r ')X2s( r ') so that it decays exponen- 
tially like Xis( r )X2s( r ) f° r large r. Thus, the second terms in the expressions (|43al) . (I43bl) for 
w X ;is(r) and r x; 2s(r) also decay exponentially, as X2s( r ) an d xL( r )' respectively. As a result, the 
self-interaction energies, r>o(ls, Is; r) and ro(2s, 2s; r), which both depend like —1/r for large r, 
dominate in the respective potentials r x; i s (r) and v x - t 2 S (r) in the asymptotic region Soo. 

As it is seen in Fig. (2b) for the Be atom, the quantities v x - n i(r)x n i(r) = F x . ni (r) + D nl . n ix n i(r) 
and v x (r)xni( r ), whose difference yields W^{ md (r), Eq. (|38l) . lie close to each other for all r. 
However, it is not straightforward to define a direct scale that could serve to estimate how small 
the potential difference v x - n i(r) — v x (r), or rather, the term W^ ,md (r) should be to make the OS 
84>nim small. Indeed, it is the ratio of the overlap integrals D n q. n i = — J °° dr'x n ii{r')W^ mA {r') 
and the orbital energy differences e n n — e n i, that, in fact, determine the expansion coefficients 
Cn'kni = D n q. n i/(e n 'i - € nl ), and, consequently, the magnitude of the OS <50 nim ; cf. Eqs. dlMU, 
(|4Q|) . Since the difference e n ' t i — e n i (with given I and n' ^ n) has the smallest magnitude for 
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n' = n + 1, we could find an upper bound for the OS norm, 

INnUI < E " vJ ^ ;n , 2 = | ll ^ ;rad " | ; (44) 

which is expressed, as it would be desired, in terms of the whole norm of W^ md (r). However, 
this bound gives values that largely exceed ||50 n ^ m || for the considered atoms; see Table HI Thus, 
it seems that, ultimately, the only fully adequate measure (in the present context) of the smallness 
of W n i" r (r) is the smallness of the norms ||50 n z m || that are generated by W^ mA {r). 



C. Properties of Hartree-Fock orbital exchange potentials 

Since the exchange-only KS orbitals a (r) = 0° EP (r) found with the exact exchange potential 
t>° EP (r) are very close to <^ F (r), the terms F x . n i(r), v x - n i(r), and D nltn i{v x ] obtained with {</>° EP } 
are virtually indistinguishable from the respective quantities F^(r), v ™^(r), D^f nl [v x ] calculated 
with the HF orbitals {0^ F } (it is true for any v x used as the argument of D n \,ni and D™ nl ). Thus, 
the combinations of these terms 

^t ra >x, {0(0 = Flf nl {r) + D™Jv x ] X ™(r) - v x (r) X ™(r) (45) 

are very close to W nl ' [v x , {0° BP }](V)- As a result, they are small for v x = v x EF (since the 
quantities W / ^' rad [t>° EP , {0° EP }] have been found to be small), and, also, by continuity, for any 
approximate potential v x close to v X EF . Therefore, basing on the numerically established prox- 
imity of the KS-OEP and HF occupied orbitals 0„/ m (r) in closed-Z-shell atoms, we conclude that 
there exists a non-empty class Vo of approximate exchange potentials v x that yield small terms 
Wni mA [ v ^ {0a F }]- m addition, we can assume that these potentials have correct, — 1/r, depen- 
dence at large r and lead to D™ H [v x ] « (since these two conditions are fulfilled by w° EP ). 

The class Vo is constituted, in fact, by all potentials v x (with correct asymptotics) for each of 
which it is possible to find constants C n \ that make terms 

Uni(r) ee *S(r) + C nlX ^(r) - v x (r) X ™(r) (46) 

small for all r and every occupied orbital x^f(r); additionally, we set Ch = 0. Indeed, this 
definition (1461) allows us to write (cf. Eq. (fT2l) ) 

POO 

D %niK\ = (Ox - v x [{0f }] ICJ = / dr X ^(r) [v x (r) X «f(r) - (r)] = 

Jo 

POO 

Cni- / dr X ™(r)U nl (r), (47) 
Jo 
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and, consequently, to express W^' iad , Eq. (|45l) . as a linear functional of U n i, namely 

POD 

W^ d K {0f }](r) = C/n/W - j£f(r) / dr'xSfCOt^CrO • (48) 



Thus, the terms ^ ;rad K, {<f>™}\ are small for any potential v x that gives small U n i, and we also 
get D™ H [v x ] « 0, due to C# = 0, from Eq. (1471) . This means that such a potential t> x belongs to 
Vo- Obviously, the appropriate constants C n i = C n i[v x ] that yield small U n i(r) for v x G Vo are not 
strictly (and, thus, not uniquely) defined with this requirement. However, according to Eq. (l47l) . 
satisfactory values of C n \ = C n i[v x ] are close to D™. nl [v x ], i.e., 

Cm[v x ] * D™Jv x ) . (49) 

Note that the small, exponentially decaying, values of U n i(r) are obtained in the asymptotic region 
for any non-diverging potentials v x , especially for those with the required, — 1/r, dependence for 
large r. 

Each approximate exchange potential v x E Vq leads to the KS orbitals <p a = a [f x ] (cf. 



Ref. 1330) that are almost identical to the HF orbitals (j)™. This can be shown by apply- 
ing the perturbation-theory argument, presented in Sec. HH to the HF equation. The orbital 
differences A<^ F = 4> a {v x } — are approximated by the first-order corrections 50^ F (r) = 
r l ^Xnf( r )^m(d, f) (where a = (nlm)) which are given by the equations (fT3~l) . ([141) where the 
KS orbitals a and energies e a are replaced with <p™ and e" F , respectively, while the perturbation 
Ah HF is used instead of 5h s . This perturbation is given by the difference h s — h HF of the one-body 
Hamiltonians entering the KS and HF equations, Eqs. ©, ©, correspondingly, so that it is the 
negative of the perturbation Ah s considered in Sec. Ill Bt Presently, we write A/ihf = —Ah s in 
the following (selfconsistent) form 

Ah w = v x -v F x [{^ F }]-Av es (50) 

and we note that the term At> es = f es [^tof] — ^esfatot] is linear in Acf)™ ps 50™ in the leading 
order. As a result, the equation (fT3l) for 5(j)^ F leads to a set of non-homogenous linear integral 
equations for the corrections 5xfJ to the HF occupied orbitals (of both spins). In these equations, 
the inhomogeneous terms (the right-hand sides) depend linearly on W^ vaA [v x ^ {0„ F }], through the 
matrix elements 

poo 

D^,Jv x } ee <0* F > X - £ x F [{0f }]|0« F m ) = - / drx^W^K {OW , (51) 

16 



o 



n' 7^ n. By solving the set of equations for 5xfJ, we can find the radial orbital differences 
Xnzfax] — Xnf ~ $Xnh which are small when all terms W„^ ad [t> x , {0a F }] are sufficiently small. 
Now, it can be claimed again (cf. Sec. III.A) that, formally, it is the norms ||50^f|| that are the 
adequate measure of smallness of W^' md [v x , {^p 7 }]. Then, the class V = V (t?) can be defined 
more precisely with the condition ||0 a [t> x ] — 0^ F || < r\ (for a = 1, . . . , N) where 77 1. 
The total energy 

E [v x ] = (^[v x }\HMv x ]) = E [{(f> a [v*)}) , (52) 



where ^[v x ] is the Slater determinant constructed of {0 a [w x ]} (cf. Ref. l33f) . is very close to 
E HF = E [{0™}] f° r an y v x G Hi ^ ue *° me orm tal proximity, (j) a [v x ] ~ <p™ . As a result, 
the energies -E[t> x ], v x E Vo, are also very close to _E[t>° EP ] since the potential v x EP minimizes 
the functional E[v x ] > E H f- The obtained relation E[v x ) ~ -E[t>° EP ] implies, by the continuity 



of the functional E[v x ] (cf. Ref. |33_|), that the potentials v x belonging to Vo are close to v x EP and, 
consequently, they are all close to each other. Simultaneously, this argument explains in a plausible 
way why the exact exchange potential v x EP itself belongs to the class Vo and, in consequence, it 
gives the KS orbitals a very close to 0^ F . 

Low magnitude of U n i(r) obtained for a potential v x E Vo implies that, within each occupied 
shell S n , the shifted HF potentials 

tiga(r)=^(r) + Ga (53) 

(I e C n = {0, . . . , 41}) lie very close to v x (r), 

?»«^), lEC n ,rES n , (54) 

and, as a result, they almost coincide with each other, 

SW«C(r), lJ'EC n ,rES n . (55) 

Similar proximity holds for the OEP potentials v x - n i(r), since they are all very close to v x EP (r) 
within their respective shells S n ; see Sec. IIII Al and Figs. [HIU 

Let us note here that since any two different exchange potentials, v x (r) and v' x (r), from the 
class Vo are close to each other, the respective constants, C n i = C n i[v x ] and C n i = C n i[v' x ], that 
lead to small terms U n i(r), Eq. (|46l) , are also close to each other. This is so because the equation 
(|54|) is satisfied for both potentials v x (r) and v' x (r), as well as for each (nl); the same conclusion 
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is reached by noting that, with Eq. (1491 ), we obtain the expression 

i>oo 

C nl [v x ] - C nl [v' x ] « D™[v x ] - D™[v' x ] = / dr [v x (r) - <(r)] [ X ^(r)] 2 . (56) 

./o 

which is small for t> x fa v ' x . In particular, by taking v' x = t> OEP we find 

C nl [v x ] fa C>*° EP ] - ^S F [- X ° EP ] - A^H (57) 

(for t> x G Vo) where the relation (l49l and the orbital proximity, a [t> OEP ] « are also applied. 
This means that the quantities D n i[v® EP ], found with the exact exchange potential t> OEP , can be 
used as the constants C„,i = C n i[v x ] suitable for all v x G Vo- Another possible set of constants 
{C n i}, which can be determined easier than D n i[v® EP ], is given by the quantities {C* LI ~ HF } = 
D^f[vf u ' HF ] found, in a selfconsistent way, for the KLI potential ti^ LI " HF = u^ LI [{0^ F }] obtained 
with HF orbitals; see Sec. MI D 31 UlI D 41 below. The two sets of constants, listed in Table HH are 
indeed very close to each other. 

A generalization of Eq. (|55l) is found when, in the expression (1461) for U n i(r), the potential 
v x (r) is replaced by #J^(r) for r G S n > according to Eq. (l54l) . the definition (l53l) is used, and the 
smallness of U n i(r) for v x G Vo is accounted for. The generalized relation reads 

+ C mX™(r) « &(r) + CW) x^ F (r) , r G £> n / (58) 

and it is satisfied for suitable set of constants {C n i} and for all indices (nl), (n'V) corresponding 
to the occupied HF orbitals, as well as for an appropriately chosen set of the shell border points 
r n,n+i- The relation (|58l is an intrinsic property of the HF orbitals (and the Fock operator), since it 
is not implied by the DFT or the definition of v 0EP , though it has been revealed here by inspecting 
the KS results for v x = t> OEP . Obviously, the relation (|58l can be rewritten as 

(00 + C ni) Xfi{r) « (^(r) + C n>v ) x« E (r) , r G S n > (59) 

so that by dividing its both sides by Xnf( r ) f° r n ' = n ' we recover the approximate equality (l55l) 
of the shifted HF orbital exchange potentials within the shell S n . 

The potentials u^(r) obtained with the occupied HF orbitals Xnf( r ) from Eqs. (1341) . (l35l) do 
not diverge for r — > oo, unlike the KS potentials f X ;n/( r )> exce pt for the HOMO one (as well as 
v x; is(r) in the Be atom). It is due to the form of the large-r dependence Xnf( r ) ~ r anl e~^ Hr with a 
common (for all nZ) coefficient /3# = a/— 2e^ F of the exponential decay and with the appropriate 
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values of the orbital- specific constants a n i, cf. Refs. |35[ [3g, |37L |38|. As a result, the following 



asymptotic dependence for the HF exchange potentials is obtained 

«£(r) = - (4 F - <%) + ^ + o f i ) (60) 



r 



where the constants q n i can differ from -1; see Appendix. This dependence is confirmed by the 
numerical results obtained for the HF orbitals found by solving the HF equation with the highly 
accurate pseudo spectral method [34]; see Fig. [51 The constant term in Eq. (l60l) vanishes only for 
the HOMO potential v™ H (r) and, in this case, we also find qu = —1. Thus, the exchange potential 
v^ F H (r) has the — 1/r dependence for large r. In consequence, it is close to the potentials v x E Vo 
not only within the region S n of the shell to which the HOMO belongs, but also in the asymptotic 
region where these potentials decay like —1/r. The asymptotic dependence 

1 fl 



< i H (r) = — + o(- , re^, (61) 
r \r J 

complements the relations (1551) . (1581) . valid within the occupied shells S n . Note that the shifted 
potential v^j(r), entering Eq. (1541) . is equal to v^ F H (r) since we set Ch = as in the definition of 
the class Vo. 

D. Accurate representations of exact exchange potential with HF orbital exchange potentials 

It has been shown above that the proximity of the individual HF and exchange-only KS-OEP 
occupied orbitals implies the relations ((551) . ((58l) satisfied by the HF orbitals. Interestingly, the 
converse is also true. Namely, assuming that the relation (1581) holds (then, the relation ((551) is also 
true) and the constants C n i which satisfy it are known, we can effectively construct local exchange 
potentials r x (r) that belong to the class Vo, i.e., which lead to small terms U n i(r), have correct 
(—1/r) asymptotic behaviour, and, in consequence, give the KS orbitals close to the HF ones. As 
it is argued above, such potentials should represent the exact exchange potential t>° EP (r) with high 
accuracy. 

1. Shell-resolved piecewise exchange potentials 

If the relations (1551 ). (1581) are fulfilled for a given set of the constants C n i, the straightforward 
way to build a potential v x E Vo is to set it equal to one of the (almost coinciding) potentials 
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u^(r), Eq. (1531) . in each occupied atomic shell S^; then, the resulting potential v x satisfies the 
relation (|54l) (which has to hold for any v x G Vo). In particular, we can choose the s-orbital 
(I = 0) potentials v^ (r) for r G 5 n , n = 1, . . . , n occ . However, within the outmost occupied 
shell S n , n = n occ , it is better to use the HOMO exchange potential v™ H (r) since it can represent 
the constructed v x (r) not only for r G S n , n = n occ , but also in the asymptotic region where it 
has the — 1/r decay (required for v x G Vo), cf. Eq. (|6TI) . In this way, a piecewise (pw) exchange 
potential 

n occ -l 

«T°(r) = £ CMCM + C(0t&M , (62) 

n=l 

is obtained; here the step-like functions #^ F '°(V) are equal to 

0(T - Oeft - r) (n < n occ ) , (63a) 
% - <n°-i) (n = riocc) • (63b) 

This construction is restricted to the case when the HOMO belongs to the outmost occupied shell, 
which is true for the closed-/-shell atoms. 

To make the potential v P w -°(r) continuous, the shell borders r n n+1 are set at the points r^ +1 , 
n = 1, . . . , n occ — 1, where its constituent potentials from the neighboring shells, S n and S n+ i, 
match, i.e., the condition 

^x H noM = ^x H n+i,o(0 (ji < n occ - 2) , (64a) 
vZo(r) = vf H (r) (n = n occ - 1) (64b) 

is satisfied for r = r™^ +1 ; we also define r™ = 0. The outer border r n „ +1 of the outmost 
occupied shell S n , n = n occ , does not have to be defined since it is not used in Eqs. (I62H63I ). 
However, if the point r n n+1 , n = n occ , needs to be determined (e.g., when we want to specify 
the region S n where the relations (1551) . (1581) or (|54l) are fulfilled for n = n occ ) it can be plausibly 
defined as the smallest of the classical turning points rJJ for electrons from the n occ -th shell; in the 
HF case, each point rJJ can be found from the condition i>^(r) + 1(1 + l)/(2r 2 ) = e^f ; cf. Eq. 
©. 

The relation (g5J (with I' = 1^2 for n' = n occ (HOMO), and V = for n' < n occ - 1 ) and 
Eq. (16Tb immediately imply that the constructed potential v x (r) = uP w '°(r) yields small U n i(r), 
Eq. (146)) . within each occupied shell S n > and also for r6S ra . Thus, the potential v P w '°(r) belongs 
to the class V and, in consequence, it is close to t> OEP (r); this conclusion is supported by the 
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numerical results plotted in Fig. Oa). Such numerical confirmation also implies that the points 



r 



„HF,0 

n,n+ 



1 can be indeed be chosen for use as the shell borders r n n+ i in the relations (|58l) . (|6"TT) . 



Another representation of t>° EP (r) is obtained by constructing a continuous piecewise potential 



?1occ 



< w (0 = E^ F W^ n) (0 (65) 



71=1 



formed from the HF shell exchange potentials 



each applied in its shell region S n . The points r™ +1 defining the shell borders are now the solu- 
tions of the continuity equation 

v^\r)=vi n+1 \r) (67) 
for n = 1, 2, . . . , n occ — 1; = 0. We denote 

p» F (r) = (2l + l)[ X ^(r)] 2 , (68a) 
Pf(r) = £ <£T(r) , (68b) 

P HF (0 = 5>f(r), (68c) 

n=l 

and the functions 9™{r) are defined like 6^ F '°(r) with Eq. (1631) where the radii r^+i 316 replaced 
by Each shell potential Ux n ^(r) is very close to the almost coinciding potentials v™ nl {r), 

I G C n , for r G S'n, Eq. (1551) . Thus, the potential t> x ™ ^ can be substituted for v^ n , in Eq. (l58l) . 
which leads to the relation 

*S(r) + C^Cr) « «<»'>(r)x*f(r) , (69) 

valid for r G 5 n / and n' = 1, . . . , n occ . It means that the terms U n i(r), Eq. (|46|) . are small for 
the potential u x (r) = w pw (r) within each occupied shell S n /. This implies that the potential y£ w 
belongs to the class Vo and, hence, it is close to t>° EP ; cf. Fig. Ob). To make this argument 
complete we note that the potential t>P w (r) is also close to v^(r) in the asymptotic region 
and, thus it has the correct, — 1/r, dependence for r G Soo (which is a property requested for 
potentials v x G Vo). Indeed, for large r, the factor p%f(r)/ p^ F (r), present in Eq. (l66l) . goes to 1 
for nl = H and it vanishes like r~ q , q > 4, for other HF occupied orbitals x^f ! ( see Eqs. (IA.1I) . 



(IA.2I) in Appendix; cf. Refs. |35|, 3jl 37, 3s[). The presented construction of t> pw (r) is restricted to 



the case of the closed-Z-shell atoms where the HOMO belongs to the outmost occupied shell. 
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2. Shell-dependent slope of the DFT exchange potential 



The slope of the exact exchange potential v ° EP ( r ) changes rather abruptly (here disregarding 
small intershell bumps) when we move through an atom, from one atomic shell to the next one; cf. 
Fig. OSIH This property can be explained by the fact that the potential t>° EP (r) is represented with 
high accuracy, within each occupied shell S n , by the potentials v™ nl (r), I G C n , and, in particular, 
by the s-orbital exchange potentials v^ F n0 (r) which exist for each occupied shell (n = 1, . . . , n occ ). 
Indeed, the slope dv^^/dr = dv^ /dr found within the shell S n for the potential v^ (r) = 
v wo( r ) + obtained with the Eqs. (1341 ), (|35l) (where the orbitals {xni} are replaced by {x^f }), 
is distinctively different from the slopes of other potentials v^, (r) within their respective shells 
S n >. It is related to the fact that the orbitals v^ ¥ nl (r) (e.g., for I = 0) corresponding to different 
atomic shells are localized at different distances from the nucleus. 

The above general argument readily applies to the Be atom. In this case, the potential 

= vf(ls, Is; r) + K 2Sils (r)v™(2s, Is; r) (70a) 

(cf. Eq. (|43al )) is very well represented, for r G Si, by the first term v (ls, ls;r); see Fig. 
I7J The other term in Eq. (|70a|) is much smaller due to the combined effect of the small ratio 

^2s,is(r) = X2s F ( r )/x?f( r ) ( we find |ft2s,is(?")| < 0.2 for < r < 0.74a.u.) and low magnitude 
of Vq F (2s, ls;r) = dr'(l/r — 1 / r')xu '( r ')x^I( r ') ( m comparison to v^ F (ls, ls;r)), which 
decreases with increasing r. We also find that, in the expression 

v™2s(r) = v% F (2s, 2s; r) + /^, 2s (r> HF (ls, 2s; r) (70b) 

(cf. Eq. (|43b| )), the term t>Q F (2s, 2s; r) clearly dominates within the shell S 2 where both the term 
t>Q F (ls, 2s; r) and the ratio K, ls>2s (r) = X^( r )/x 2 I( r ) decay exponentially for the Be atom; cf. 
Fig. |7j Thus, we obtain the relation 

dv wi( r ) „ dv™(nl,nl;r) Q ni (r) 
dr dr r z 

which holds for both nl = Is and 2s; here Q n i(r) = J Q r dr'[x^(r')} 2 . As a result, we conclude that 

the derivative dv™ ls {r) / dr « —3.65 a.u. at the point r = r* s = 0.37 a.u G Si where Qu(r) = 0.5 

has the magnitude (approximately) (r^/r^) 2 = 44.2 times larger than the slope dv™ 2s {r)/dr m 

—0.0826 a.u. at the point r = r| s = 2.46 a.u G S 2 where Q 2s (r) = 0.5. These estimates agree 

well with the values —3.70 a.u and —0.0846 a.u of dv° EF (r)/dr at the points r = r* s and r 

respectively; the ratio of these two slopes is 43.7. 
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3. KLI- and LHF( CEDA )-like potentials constructed from the HF orbitals 



The KLI-like potential v^ l (r) can be defined for the HF orbitals {4>™ m } and the constants 
{C n i\ by substituting them for {<p n i m } and {D aa } = {D nl . nl }, respectively, in Eqs. (Q3I2D- It 
takes the following form 



-KLI-HF 



[{C nl }} (r)^ 1 [{x^}ACm}] (r) = 

occ 

£(2Z + 1) [F™(r) + C nlX «!(r)} X » F (r) (72) 

P ^ ' nl 

where, for the closed-Z-shell atoms, the quantities Xnf an d C n i are indicated as the effective argu- 
ments of v^ 1 . This potential can also be expressed in terms of the HF orbital exchange potentials, 

-x™ [{C nl }] (r) = £ ^(r) . (73) 

nl P ^ ' 

It can be argued that the potential i;f LI '" HF (r), Eq. d72l . is close to ^ LI [{x° EP }, {-D°^}](r), and, 
consequently, also to v ® EP (r) (cf. Sec. Ill CI) , because the HF orbitals Xnf( r ) near ly coincide with 
the KS-OEP orbitals Xnf P ( r )> while the constants C n i satisfying the relation (1581) are very close to 
D nfja = D n i ]n i{v° EF }; cf. Eq. (|57l> . However, the high-quality of the KLI-like potential i^ LI " HF (r) 
is, in fact, a direct consequence of the relation (1581) revealed for the HF orbitals. Indeed, this 
relation immediately implies that the potential given by Eq. (1721) is close to v™ n , v (r), V E £ n >, for 
r E S n >,n' — 1, . . . , n occ . This means, in particular, that the potential tj| aj ~ HF (r) is close to t>P w '°(r) 
within each occupied shell S n i so that it also yields small terms U n i(r) there (for any (nl) E occ). 
For large r, the potential #f LI " HF (r), given by Eq. (1731) (with C# = 0), becomes close to vf F H (r) 
so that it decays like —1/r (see the discussion for t>P w (r) above). These properties of the potential 



v. 



KLI-HF ; 



imply that it belongs to the class Vo and, in consequence, it is close to v ° EP , cf. Fig. (6fc). 



In particular, this is true for the KLI potential 

^KLI-HF^) = ^KLI-HF (r) = ,KLI ^j] (r) (?4) 

calculated with Eq. (1231 for the HF orbitals and the constants C*p I_HF that are found from their 
self-consistency condition 

D^Jv™-™} = CT HF (75) 

I — I 

given in Ref. [5]. To show this, let us express them as the sum 

CT™ = C nl + AC nl , (76) 
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AC nl - £(2J' + 1) ^ dr ^^^ AC nq> = - dr X ™(r)U nl [v™-™}(r) (77) 

„/;/ JO P V) JO 



-jKLI-HF 



where the constants C n \ satisfy the relation (l58l) . Then, we obtain, from Eq. (l47l) (the first line) 
and Eqs. (|72|) . (1751) . the following set of linear equations for AC„i 

[xST(r)^(r)] ! 

n'i' 

(n = 1, . . . ,n OC c, I E C n ; nl ^ H) where the right-hand side includes the potential 
-KLi-HFj-|^ jj^ gq (f72l) . The set of equations (1771) remains satisfied when a common constant is 
added to each AC n j. Therefore, to make this set well-defined, we put ACh = and, simultane- 
ously, exclude the equation for nl = H from the set (then, we find C^ LI " HF = C# = 0). Since the 
terms U n \ [■u x CLI ~ HF ] (calculated for C n \ satisfying Eq. (1581) ) are small, the corrections AC n i obtained 
by solving the equations (fTTl) are also small. This means that the potential v^ 1 '™ is very close to 
-kli-hf^^ |j e y Q an( ^ m conS equence, this potential itself belongs to the class V . 

The KLI condition (1751) can also be satisfied by minimizing, with respect to the constants {C n i}, 
the function 

occ occ 2 

g ({C nl }) = J2 \MT - «W - GjO 2 = E( 2/ + !) / dr (^N(r)) (78) 

where we put v x = vf L1 ' HF [{C n i}], Eq. (l72l) . and a = (nlm); a similar expression leads, after 
minimization, to the selfconsistent constants {D a b} for the LHF (CEDA) approximate potential 



L3 



4(j]. To avoid the presence of an arbitrary common constant that can be added to all C n \ = 
" HF ( smc e such addition does not change the value of g ({C n i})), we again set C^ LI ~ HF = in 
Eq. (1781) . The function g ({C n i}) attains very small value for the constants {C n i} that satisfy the 
relation (1581) since they lead to small terms £/nj[?; x ] for t> x = , u^ LI " HF [{C n j}]. The set of constants 
{C n i} = {CnZ LI ~ HF } that minimizes the function (1781) have to yield even lower value of g ({C n i}), 
and, in consequence, they should also give small terms U n i [v x ajI " HF ] . Thus, we can conclude again 
that the corresponding potential v^ 1 '™ belongs to Vo- 

By extending the arguments presented above for the KLI potential we can show that the high 
accuracy of the LHF (CEDA) approximation is also directly explained by the revealed properties 
of the HF orbital exchange potentials. Let us first note that the potential vf L1 [{xjjf }, {C?a}] is a 
special case of the LHF-like potential 

^x LHF - HF [{{Cn>i,m}\ (r) = £ X LHF M F }, {C nH , nl }] (r) = 
(21 + 1) [<WS(0 + Cn>i,mX^i{r)] j£f(r) 

(79) 
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occ 



p HF (r) ^— ' 

r v ' nLn'l 



calculated with Eq. (l29l) for the HF orbitals and the constants C n 'i >n i = 5 n / n C n i. We can now solve 



the LHF self-consistency condition 1117 



1H 

n [ LHF-HF] _ ^yLHF-HF /o m 

where 



- x LHF - HF (r) - ^ X LHF 



[{)8h{<%gF}] (0, (81) 

by expressing C^,^ HF as 5 n ' n C n ; + AC^/^j. Then, the corrections AC n i^ n i satisfy a set of linear 
algebraic equations (similar to Eq. (1771) ) where the right-hand sides are given by the integrals 
— J °° drXnn( r )U n i(r); we also set ACh,h = 0. The terms U n i(r) are small since they are cal- 
culated here for v x = v^ 1 '™ [{C n i}} obtained with the constants C n i satisfying the relation (l58l) . 
Thus, the resulting corrections AC n >i >n i are also small. This implies that the potential t> LHF " HF (r) is 
close to i;^ 1 "™ [{C n ;}] (r) and, as a result, it also gives small U n i(r). In effect, the LHF exchange 
potential v LHF " HF , obtained with the HF orbitals, belongs to the class Vo and it is close to v ° EP . 



4. Comparison of different approximate representations of exact exchange potential 

The constants {C*P" HF } obtained with the KLI selfconsistent condition (1751) have been shown 
to differ only by small corrections {AC n i} from from any set of constants {C n {\ satisfying the 
relation ((58l) . This property combined with Eq. (1571) explains the small magnitudes of C* LI ~ HF — 
D n i;ni[v® EP ], cf. Table|IIl It also implies that the constants {C*^ I ~ HF } themselves satisfy the relation 
(l58l) so that they can indeed be used in construction of the approximate potentials discussed in Sec. 
IIII Dl In particular, as it is already mentioned above (see Fig. [6]), the potentials uP w,0 (r), vP w (r), 
built entirely with the HF orbitals and the constants C n i = C^ LI " HF , are found to be very accurate 
representations of the exact exchange potential t>° EP (r). The obtained quality of its approximation 
is almost the same as for the potentials w^ LI " HF (r), Eq. (1741) . and ujP^r), Eq. (fT9l) , the latter of 
which is built of the KS-OEP orbitals and it is the dominant part of v° EP (r), Eq. (TT"8~1) . However, any 
of the four approximate potentials fails to reproduce well the characteristic bumps of v ° EP (r) at the 
shell borders; cf. Figs. OS HI [61 Thus, it is the minor part of the exact exchange potential, namely, 
its OS term v ° s (r), Eq. (f22]) . depending linearly on 5(J) a (r) (a = 1, . . . , N), that produces these 
local maxima of v® EP (r). This means that the intershell bumps of v° EF (r) are the consequence 
of the finite (though very small) differences o (r) — 0Q F (r) « S(p a (r) between the KS and HF 
occupied orbitals. 
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The potentials v P w '°(r) and t>P w (r) are expressed, in each atomic shell S n , in terms of the orbital 
exchange potentials v™ nl (r), I G C n , that correspond to this shell only. This feature makes these 
two representations of the exact exchange potential be significantly different from the KLI-like 
potential, Eq. (|73l ). Indeed, the latter depends, within each shell S n , on all potentials v™ nll (r), 
corresponding to both the same (n' = n) and other (n' ^ n) shells. In consequence, the KLI-like 
potential Eq. (1731) . rewritten as follows, 

^occ HF / \ 
n'=l ^ ^ ' 

is given for r e S n not only by the respective shell potential Vx (r), but it also expressed there 
by the potentials v x (r) which correspond to other shells (n ^ n) and can be calculated for any 
r with Eq. (|66l) . The fact that, despite its significantly different structure, the potential w^ LI " HF (r) 
(calculated with appropriate constants {C n i}, e.g., {C^ LI ~ HF }) is very close to wP w '°(r) and v^ w (r), 
and ultimately also to u° EP (r), has been shown above to result from the relation (1581) which holds 
for all nl E {occ} within the occupied shells; for large r the three approximate potentials and the 
exact one nearly coincide with each other due to Eq. (|6TT) . 

Finally, let us note that the potentials vf'°(r) and v£ w (r) are identical for the Be atom since, in 
this case, there are only two occupied orbitals, one in each of the two shells. 

E. Energy shifts. Step-like structure in the response part of exchange potential 

It is known that the energies e^f of the KS-OEP occupied atomic orbitals are higher than the 
corresponding HF energies (except for the HOMO energies which are nearly equal in the two 
schemes); see Table HB The differences Ae n i = e^f — are non-negative and, for given I, they 
are the larger the lower shell index n is. The present results shed some light on these numerical 
findings as it is shown below. 

Since the KS-OEP shifted orbital exchange potentials v x - n i(r) and v X ;n+i,i{ r ) ( as wei l as me 
respective HF potentials) match quite closely at r = r n n+1 (cf. Figs. mSIU), we find 

AD n i in i = D n i in i — D n+ x,l;n+l,l ~ ^x;n+l,«( r ",n+l) — V x;nl (r n ,n+l ) > . (83) 

The latter inequality results from the mathematical structure of the Fock operator £> F (r), Eqs. 
(|33~T) . (I34al) . (I34bl) presumably because the orbital Xni( r ) is localized closer to the nucleus than 
Xn+i,i(r)- This argument is certainly valid for the Be atom. In this case, the terms proportional to 
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t>o(ls, 2s; r), which are present in Eqs. (I43al) . (I43bl) . are negligible at the point r = r i2 = = 
0.954 a.u. where X2s( r )/xis( r ) = 0.69 (see Fig. [7]); as a result, we have 

Wx;2s(ri 2 ) - v x - ls (r 12 ) « w (2s, 2s; r i2 ) - u (ls, Is; r i2 ) . (84) 

The latter difference can be found by integrating the equation (TTTT ) (here for the terms v (nl, nl; r) 
defined with the KS orbitals) , 

rri2 i 

t;o(2s,2s;r 12 )-^(ls,ls;r 12 ) = / dr'— - (Q ls (r') - Q 2s (r')) > , (85) 

./oo v / 

and it is positive since the relation Qi s {r) > Q 2s (r) holds for any r. Let us note here that the 
approximate relation (f8~3l) is not satisfied very tightly for the closed-Z-shell atoms other than Be 
since the differences Af x;n+ i i /(r) = t> X ;n+i,/( r ) — v x;n i(r) change quite rapidly around r = r nyTl+ i 
(due to very different slopes of the orbital exchange potentials from the neighboring shells; see 
Figs. Ute), 121 SI) while the point r where the potentials v x;n i{r) and w X ;n+i,/( r ) intersect slightly 
differs (except for the Be atom) from the shell border r n n+1 = r^ F n+1 (defined in Sec. IIIID II) : cf. 
Fig. H£d). However, as it is seen in Table HH the differences Av x;n+ i t i(r n >n+ i) have quite similar 
value and definitely the same sign as the corresponding constants AD^. 
Further, we can express D n i tTl i as follows 

h(l)-l 

D n l t nl = D 'n(l)l,n(l)l + AD n 'i >n 'i (86) 

n'=n 

for n < n(l) where the symbol n(l) denotes the largest shell index n among the KS-OEP occupied 
orbitals Xni(r) with given orbital number I. Thus, according to Eq. (l86l) and the inequality (1831) . 
the energy shift Ae n i w D n i )Tl i grows with decreasing n and, consequently, it is positive for n < 
h(l) provided the shift D^/niMhi is non-negative. The latter condition is satisfied by the HOMO 
shift D HH which vanishes. For other orbitals Xau)h the relation Dfb(i)iMi)i > is established 
numerically but understanding its origin needs further study. 

The revealed representation of the exact exchange potential t>° EP (r) with the (both HF and 
KS) orbital (or shell) exchange potentials does not result from the characteristic properties of its 
response part 

v ™v(r) = v° EP (r) - vf(r) = vf(r) + v° s (r) . (87) 



This term has been found numerically Il27ll to have a nearly step-like dependence on r where each 



step corresponds to an atomic shell. The main part of v x esp (r) is the energy-shift (ES) term 

ES / \ _ EZ C ( 21 + VDnl,nlX 2 nl(r) 

x[) £^ + l)x^(r) {m 
27 



obtained from Eq. (TfTI) . The step-like r-dependence of t>£ esp (r) fa t> ES (r) is briefly explained in 



Ref. 
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by noting that within a given shell S n the orbitals Xn'i'(r), n' ^ n, corresponding to other 
shells, are small so that they can be neglected in Eq. (f88T ). This argument can be supplemented 
by the numerical fact that the different occupied orbitals Xni( r ) (I G £ n ) from the n-th shell have 
similar shapes and magnitudes within the respective shell region S n . 



IV. CONCLUSIONS 



In summary, we find that when, for each HF orbital, a suitably chosen (orbital-specific) constant 
shift is added to the Fock exchange operator in the HF equation, the electrons occupying different 
HF orbitals are subject to very similar local exchange potentials (as well as the total ones) within 
the atomic regions where the radial probability densities of the respective orbitals are substantial. 
This proximity is particularly tight for the shifted exchange potentials of the orbitals that belong 
to the same shell and it holds in the region of this shell. Thus, the occupied HF orbitals are only 
very slightly disturbed when the orbital- specific shifted exchange potentials are replaced in the HF 
equation with a common exchange potential that lies very close to them within their respective 
shell regions; simultaneously, the corresponding orbital energies change considerably since the 
applied shifts are quite sizeable. As a result, the DFT exact exchange potential t>° EP (r) (obtained 
in the OEP approach by minimizing the HF-like total energy expressed in terms of the KS orbitals 
coming from a common local total potential) is very well represented in each shell with the HF 
shifted orbital exchange potentials from this shell, and, even slightly better, with their weighted 
average - the shell exchange potential, Eq. (l66l) . Thus, the shape of the DFT exchange potential in 
atoms, as well as its strongly shell-dependent slope, are, in fact, determined by the r-dependences 
of the individual HF orbital exchange potentials within their corresponding shells. 

The revealed properties of the shifted orbital exchange potentials result from the more general 
relation (|58l ) satisfied by the Fock exchange operator and the HF orbitals. Thus, it is in fact this re- 
lation that explains the outstanding proximity of the HF and KS orbitals in the closed-/-shell atoms 
as well as the high-quality of the KLI and LHF(CEDA) approximations to the exact exchange 
potential t>° EP . However, since these approximations are expressed in terms of the exchange po- 
tentials of all occupied orbitals (at a given point r), one of qualitatively new achievements of the 
present work is showing that the potential v° EP (r) can be represented, with equally high accuracy, 
by the (HF or KS) individual shifted orbital exchange potentials within their corresponding shells. 
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An intermediate stage between these two types of representation is obtained with the piecewise 
function formed with the shell exchange potentials. It is also shown that the positive values of 
the differences e^f — between the energies of the respective KS and HF orbitals, as well as 
their increase with decreasing n are related to the differences between the orbital exchange poten- 
tials from neighboring shells at the shell borders. Finally, it should be stressed that the presently 
obtained shell-resolved mapping between the HF orbital exchange potentials and the DFT exact 
exchange potential is not related to the previously established step-like structure of the response 
part of the exchange potential. 



APPENDIX: ASYMPTOTIC DEPENDENCE OF HARTREE-FOCK EXCHANGE ORBITAL 
POTENTIALS 



In the atomic region outside the occupied shells, the HF orbitals have the following asymptotic 



dependence Q35 



36, 



37, 



38] 



-far 



(A.l) 



where the coefficient (3h = a/— 2e^ F is common for all nl while the constants b n i and the powers 
a n i are orbital- specific. The largest a n i is found for the HOMO and it is equal to an = 1/ Ph for 
neutral atoms. For other HF orbitals the powers a n \ depend on the orbital number I, i.e., 



«ni = OLH — 3 

®nl — OLH - 2(/ m in + 1) 

®nl = OLH — \l - Ih\ - 1 



(/ = WO, nl^H), 
(l = l H = 0, nl ^ H), 



(A.2a) 
(A.2b) 
(A.2c) 



Here, In denotes the HOMO orbital number and / m i n is the smallest non-zero I within the set of 
the occupied HF orbitals in a given atom. The above asymptotic dependence (IA.1I ) is valid for all 
atoms other than Be. 

In the asymptotic region, the HF hamiltonian /i HF (r), Eq. ©, is dominated by the kinetic and 
the exchange terms since, for a neutral atom, the sum v ex t(r) + f^ F (r) decays exponentially (as 
e -2fary f or laj-gg r Thus, the HF radial equation has following asymptotic form 



2r 2 



(large r) 



(A.3) 
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and, by dividing its both sides with % 





1(1 + 1) 



2r 2 




) 



(A.4) 



When the asymptotic dependence (IA.1I) of the orbital xfa( r ) 1S applied, the general asymptotic 
form (l60l) of the HF exchange orbital potentials is found. 

By using the explicit expression for u^j(r) (given by Eqs. (I34al) . (I34bl) . (T331) with the HF or- 
bitals), one readily finds the asymptotically dominating term (—l/r)v Q (H, H;r = oo) = — 1/r in 
the HOMO exchange potential vf F H (r); cf. Eq. ( |6T| ). The same term, (— l/r)v (nl, nl; r = oo) = 
—1/r, is present in the asymptotic dependence of any potential v™ nl (r), but, for nl ^ H, it also 
includes other terms which are proportional to 1 jr or tend to constant values for r — > oo (the latter 
contribute to the constant term — (e^ F — e^f) in Eq. (l60l) ). For instance, the potential v^ p (r) con- 
tains the terms proportional to r~ 3 X3 P (r)v 2 (3p, 2p; r)/x2 P (r) and r~ 2 X3 S (r)v 1 (3s, 2p; r)/x2 P (f) 
which depend like (c x + c 2 /r) and c 3 /r, respectively, for large r; here c l9 c 2 , and c 3 are constants; 
these asymptotic dependences can be derived using Eqs. (IA.1I ) and (IA.2I) 

For the Be atom, the two occupied s orbitals decay as Xno( r ) ~ r 1 / /3 e~ /3r (n = 1,2) where 
(3 = a/— 2e"o- Thus, according to Eq. (IA.4I) . the potentials u^ ( r ) = 1, 2) vanish at r — > oo. 
They have the same asymptotic dependence (— l/r)v (nO, nO; r = oo) = —1/r, which results 
from Eqs. (I70al) . (I70bl) and the definition of the functions v^F(n'l', nl; r), Eq. (|34bl) . 
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TABLE I: The norms of Acj) n i m = cjQ m - (\> n[m , 8<j> n i m , A0 nim - and (in the last column) the 

upper bound of ||<5</> n j m ||, Eq. (I44l ). for the occupied orbitals 4> a = 4> n i m in the Be and Ar atoms; see text 
for details. 



atom 


orbital 


||A^ TO || 


WHnlmW 


||A0 n ; m - {-5(t> n i m )\\ 


||^ ;rad ||/|6 n+li/ -6 nZ | 




(nl) 


(xlO" 3 ) 


(xl0~ 3 ) 


(Xl0- 3 ) 


(xl0~ 3 ) 


Be 


Is 


6.0890 


6.6865 


0.6253 


11.2134 




2s 


5.7655 


6.3021 


0.6416 


242.9996 


Ar 


Is 


1.2594 


1.2752 


0.0305 


3.3526 




2s 


6.2281 


6.5057 


0.2929 


22.0419 




2p 


4.3019 


4.5323 


0.2467 


82.4518 




3s 


5.8187 


6.4366 


0.8003 


122.1715 




3p 


4.3474 


4.5782 


0.3428 


242.2264 
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TABLE II: The HF and exchange-only KS (OEP) orbital energies, e n i, e^f, in the Ar atom. The difference 



KLI-HFl 



e n « - compared with the constant shifts D°f p = D n i{v° EP ] = 5e nl and {C^ 1 '™} = D^f[v^ 
obtained with the OEP and KLI-HF exchange potentials, respectively. The constants AD^ n+1 = D®f p — 
.Dn+i l compared to the differences Av^) n n+ i(r) = v x - n+ i : i(r) — v x;n i(r) of the KS-OEP orbital exchange 
potentials at the shell borders r = r^ n+1 . See text for details. Note that if the point r^J was moved by 
just 0.01 a.u., to r = 0.137a.u. for the Ar atom, the considerably modified value Avj^(r) = 2.98hartree 



would be obtained . The HF orbitals and their energies used in the calculations are taken from Ref. 



30. All 



enerj 


*ies and radii are given in hartrees. 














atom 


nl 


HF 

t nl 






d oep 

nl 


r^KLI-HF 


ad' 1 ' . . 

n,n + l 


x;n,n+l 


r HF 

' n,n + l 


Be 


Is 


-4.732 669 


-4.125 699 368 4 


0.606 969 


0.606 401 428 6 


0.562 484 


0.606 


0.607 


0.954 




2s 


-0.309 269 


-0.309 227 738 5 


0.000 041 


0.0 


0.0 








Ar 


Is 


-118.610 349 


-114.452 154 608 6 


4.158 194 


4.156 319 209 3 


4.153 224 


2.991 


4.426 


0.127 




2s 


-12.322 152 


-11.153 224 215 2 


1.168 928 


1.165 666 206 9 


1.126 130 


0.988 


1.031 


0.729 




2p 


-9.571 464 


-8.733 757 145 4 


0.837 707 


0.837 222 865 2 


0.764 760 


0.837 


0.963 


0.729 




3s 


-1.277 352 


-1.099 246 843 1 


0.178 105 


0.178 063 552 5 


0.180 419 










3p 


-0.591 016 


-0.590 751 487 8 


0.000 265 


0.0 


0.0 
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FIGURE CAPTIONS 



Fig. 1. OS norm square ||50 a || 2 (grey bars) and the contributions c\, lnl (stacked bars) to it 
from bound states (f) n >i m , for the occupied states <\) n [ m in the Ar atom; the contributions from the 
occupied states are marked with the hatch patterns; the Is bars are magnified by the factor 20. 
The results are obtained in the exchange-only KS-OEP scheme. 

Fig. 2. (a) KS-OEP radial electron density p (per spin) and (b,c) the term F x . ni + D n ^ nl x n i 
(dashed and dotted lines) compared to v x EF Xni (solid lines) in the Be atom, (nl) = Is, 2s. (d,e) 
The potentials v x EP (solid line), v x; i s (dashed-dotted line), v x . ls (dotted line), v x . 2s = v x -2 S (dashed 
line), t>^ I - HF (long-dashed line in the insert (e)). The HF radial electron density p HF and the HF 
potentials v^, v™ nl , nl = Is, 2s, follow p, v x - n i and v x - n i, correspondingly, within the figure 
resolution. The up and down arrows mark the points rf E and r™", respectively. 

Fig. 3. (a) KS-OEP radial electron density p (per spin) and (b,c) the potentials v x EP (solid 
line), v x - n i (dashed and dotted lines) in the Ar atom, (d) The differences Av x;n i = v x - n i — v x EF 
(dashed lines) and Av x - n i = v x - nl — v x EF (solid lines), each shown within the r-interval in- 
cluding the corresponding shell S n and slightly overlaping the neighboring shells (SVi-i and/or 
S n+ i). The HF radial electron density p™ and the HF potentials v™ nl as well as the differences 
Av x E nl = v™ t - v° EP and Av^, = v^ nl - v° EP follow p, v x . nl , Av x - nt , and Av x - nt , correspondingly, 
within the resolution of the respective figure. The up and down arrows mark the points +1 and 
r™ n , respectively. 

Fig. 4. Results for the Zn atom; the description of the panels (a)-(d) as in Fig. [3] The HF 
quantities p HF , v x E nl , Av x E nl , Av™ nl follow p, v x . nl , Av x - nh and Av xynb correspondingly, within the 
resolution of the respective figure. 
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Fig. 5. (a,b,c) Asymptotic dependence of the potentials v™ nl (solid lines) and v x - n i (dotted lines) 
compared with the HF asymptotic limits, equal to — (e# F — e^f) (horizontal dashed lines) , Eq. 
(|60l ), in the Ar atom, (c) The HOMO exchange potentials and v x - 3p (which follow each other 
within the figure resolution) are compared with the — 1/r (dashed line) asymptotic dependence of 
t> OEP . The results are obtained with the KS-OEP and HF orbitals calculated, with high accuracy, 



by using the pseudo spectral method DUJ, l34£] . Note that the divergence of v™ s (r) seen in the panel 
(a) results from the node of the HF orbital x?f( r ) at r = 1.09 a.u.; this node is also present in 
X H J(V) calculated with the Slater-type-orbital expansion given in Ref. 



Fig. 6. Differences between approximate and exact exchange potentials: (a) t> pw - — v° , (b) 
v r - v x E? > (c) ^- HF - v° EP , (d) vf L \r) - v° EP = -v° s (Eqs. CD§, W)', see text for details. 
The dashed lines correspond to v X EP — v X EF = 0. The up arrows mark the points r^ F +1 , n = 1,2. 



Fig. 7. (a) KS-OEP orbital exchange potentials v x - ls , v x - 2s (solid lines), Eqs. (|43a|) , ( |43b| ), com- 
pared with the contributing functions v (ls, Is; r), t> (2s, 2s; r), v (ls, 2s; r) = t> (2s, Is; r) (dot- 
ted and dashed lines), Eq. (|34bL and (b) the ratios X2s/Xis, Xis/X2s for the Be atom. The HF 
potentials t>^ F , Eqs. (I70al) . (I70bl) . functions t> HF (nO, n'0; r) and ratios Xno/Xn'o ( n > n ' = h 2) fol- 
low the corresponding KS-OEP quantities within the figure resolution. The up down arrow marks 
the difference Av^) nn+1 (r) = v x - n+ i^{r) — v x - n i(r) at r = r" F = 0.954 a.u.; it is very close to 
L>oep = 0.6064 hartree; seeTableHH 
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